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Abstract 

We calculate mesonic two-point functions in the e-regime of two-flavor QCD on the lattice with 
exact chiral symmetry. We use gauge configurations of size 16^ X 32 at a ~ 0.11 fm generated 
with dynamical overlap fermions. The sea quark mass is fixed at around 3 MeV and the valence 
quark mass is varied in the range 1-4 MeV, both of which are in the e-regime. We find a good 
consistency with the expectations from the next-to-leading order calculation in the e-expansion 
of (partially quenched) chiral perturbation theory. From a fit we obtain the pion decay constant 
F = 87.3(5.6) MeV and the chiral condensate Y.^^ = [239.8(4.0)MeV]3 up to next-to-next-to- 
leading order contributions. 
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I. INTRODUCTION 



In the low energy limit, the dynamics of Quantum Chromo dynamics (QCD) is domi- 
nated by the pion fields that appear as pseudo-Nambu-Goldstone bosons associated with 
the spontaneous breaking of chiral symmetry. Although chiral Perturbation Theory (ChPT) 
is a powerful effective theory in understanding their interactions, it has many parameters, 
the so called 'low energy constants (LECs)', which should be determined either from experi- 
mental data or preferably from direct calculations based on the underlying theory, i.e. QCD. 
The leading order of ChPT is completely controlled by the two LECs, the chiral condensate 
S and the pion decay constant F, while at higher orders there are increasing number of 
LECs. 

Numerical simulation of lattice QCD offers the most promising approach to achieve the 
direct calculation of LECs. In fact, the recent advances in the simulation techniques allow us 
to calculate, for instance, the pion decay constant to a remarkable precision. However, such 
calculations may contain non-trivial systematic effects, since the infinite volume limit must 
be taken before the chiral limit is approached. The violation of chiral symmetry of lattice 
fermions also seriously complicates the analysis, because ChPT itself must be modified to 
incorporate the explicit breaking of chiral symmetry 

Recently, an alternative approach has been proposed, that is the lattice calculation in 
the e-regime of ChPT . In this approach, the lattice simulation is performed 

near the chiral limit at a fixed volume V. Finite volume effect becomes prominent due to 
the long distance correlation of the pion fields, which can be treated in a systematic way 
within ChPT. Of particular interest is the region where the pion correlation length, or the 
inverse pion mass exceeds the size of the box L 

^ , (1) 



Aqcd 

with Aqcd the QCD scale. In this e-regime, the zero-momentum mode of the pion fields has 
to be treated non-perturbatively, and the contribution from non-zero momentum modes is 
expanded in a new parameter e: 
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where p denotes pion momentum and Acut is a cutoff of ChPT. With this expansion, the 
volume and topological charge dependence of the chiral condensate, meson correlators, etc., 
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can be written with the same low energy constants as in the infinite volume. 

Simulating lattice QCD in the e-regime has several advantages over the conventional 
approach. First of all, the infinite volume limit is not necessary when approaching the 
chiral regime. In the p-regime, on the other hand, it is not clear at which masses one can 
safely apply the chiral expansion to extrapolate lattice data to the physical up and down 
quark masses. This is a question that depends on the quantity of interest, and therefore one 
potentially needs to work on a fairly large volume lattice to ensure the condition rrtT^L ^ 1. 
In the e-regime, pion mass is made arbitrarily small at a finite volume. The e-expansion 
requires AirFL ^ 1, which does not strongly depend on the quark mass. 

One may wonder that the computational cost to simulate arbitrarily light sea quark could 
be prohibitively high, but it is not the case because the lowest-lying quark eigenvalue stays 
finite at the order of 1/T,V for a given volume V (except for the exact zero-modes). In the 
e-regime the lowest eigenvalue is even lifted by about a factor of A'^^, the number of fiavors, 
according to the chiral random matrix theory. 

Another major advantage of the lattice calculation in the e-regime is the prominent 
dependence of physical observables on the topological charge and sea quark mass. The best 
known example is the analytic predictions for the lowest-lying eigenvalues derived from the 
chiral Random Matrix Theory. Utilizing these, precise determination of the low energy 
constants has been attempted through direct calculations of the low-lying Dirac eigenvalues 



mainly in quenched QCD 



Another possibility is to study hadron correlators in the 



e-regime, which has also been carried out by several groups in quenched QCD 
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161 . We extend these works to unquenched QCD in this paper. However, this is not 



an easy task since the exact chira! 
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symmetry is essential in the e-regime. We therefore use 
3l| which realizes the exact chiral symmetry on the lattice 
3l|. The overlap fermion is much more difficult to simulate than other lattice fermion 



;he overlap-Dirac operator 
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formulations for both algorithmic and computationa' 



The recent series of work by JLQCD collaboration 17l . ll8l.ll9 | has opened a new 



reasons. 



possibility 



of simulating unquenched QCD in the e-regime (see also [20|, l2l|, [2^, l23|, [24, l25|, l26|, [27 1 for 
other exploratory studies). We performed two-fiavor dynamical overlap fermion simulations 
with the quark mass near the chiral limit, on a 16^ x 32 lattice at a lattice spacing a ~ 0.11 fm 



(determined with tq ~ 0.49 fm 3^ as an input). We use the Iwasaki gauge action 
with extra Wilson fermions and ghosts to fix the topological charge 35 



36 



33, 



34| 



37|. The sea 
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quark mass m is around 3 MeV, which is well within the e-regime. Since the eigenvalue of 
the hermitian overlap-Dirac operator is bounded from below, numerical simulation is stable 
even with such a small quark mass. Comparing the Dirac spectrum with the predictions 
of the chiral Random Matrix Theory, we extracted the value of the chiral condensate at 
the leading order of the e-expansion as = [251 (7) (11) MeV]^ where the second error 
is an estimate of the systematic error due to the next-to-leading order (NLO) effects in the 
e-expansion. 

In this paper, we use the same set of gauge configurations in the e-regime to calculate 
the meson correlators in various channels. The analytic predictions of ChPT for the pseudo 
scalar, scalar, axial- vector, and vector channels are known to NLO in the e-expansion 
which are recently extended to the partially quenched ChPT 40|, 
predictions to extract S and F at the NLO accuracy. 

This paper is organized as follows. In Section [TTl we review the (partially quenched) ChPT 
predictions for the meson correlators. The set-up of the numerical simulations is given in 
Section [ml In Section HVl we measure the axial-vector and pseudo-scalar correlators to 
extract F and S. Then, some consistency checks are done using other channels and partially 
quenched correlators. Comparison of the result for S is also made with that from the Dirac 
eigenvalue spectrum. Our conclusions are given in Section IVII 
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39|, 



4^. We use these ChPT 



II. (PARTIALLY QUENCHED) CHIRAL PERTURBATION THEORY IN THE 
e-REGIME AT FIXED TOPOLOGY 



In this section, we briefly review the results for the meson correlators calculated within 
(partially quenched) ChPT. For the full details we refer the original papers 38|, l39|, |41| . Here 
we consider valence quarks with a mass rriy and A'^^ = 2 degenerate sea quarks with a 
mass rris, both in the e-regime. 

As a fundamental building block for the later use, let us define the partially quenched 
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44j at a fixed topological charge u. In addition to the 



zero-mode partition function 
Nf = 2 physical quarks (of mass m^), a pair of a valence quark (of mass my) and a bosonic 
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quark (of mass rrib) is introduced for the partial quenching; 

/ K^iflb) LifJ'v) Lif^s) Iu-l{f^s)/fJ's\ 

f^b^i^+2ifJ'b) f^lh+2ifJ'v) fJ''ilu+2ifJ's) f-l'slu+lif^.s) 

\ -f^lK^+sif^b) fi%+3{f^v) /i^/^+3(/i.) fi%+2{fis) J 



det 



(3) 



where Hh = rubTy , /i^ = m^'EV, and fig = rris'EV. K^s and J^'s are the modified Bessel 
functions. Note that, in the hmit /i^ — > ([3]) reduces to the zero-mode partition function 
of the full Nf = 2 theory: 



lim = Z'f{iJis) = det 



Then the partially quenched chiral condensate at finite V and u is given by 



(4) 



- lim 

-1 



lnZ^^{fib\fi'v,fJ's) 



X det 



d^AfJ'l^<^u+2M) /i^4+2(/i^,) fi'ilu+2{f^s) fisL+lif^s) 
\ -di,^{fllK^+3{fi^)) flll^+sifi,,) fl^I^+sifls) ^^llu+2{f^s) J 



(5) 



It is not difficult to see that in the fiy fig limit, the partially quenched condensate reduces 
to the one in the full theory. 



lnZr(/is). 



S S 2 

In the following, we will also use a second- derivative 

AEPQ(/i^,/i,) _ lim^,^^„9^,9^„ZPQ(/ib| 



(6) 



(7) 



First we present the two-point correlation functions of the flavored pseudo-scalar and 
scalar operators, P°-{x) = q{x)T°'^^q{x) and S°'{x) = q{x)T°'q{x), where denotes a gener- 
ator of SU{Ny) group that the valence quark field q{x) belongs to. 
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The expressions for the correlators (P°'(x)P°'(0)) and {S°'{x)S°'{0)) in the partially 
quenched ChPT is known to 0{e^) (no sum over a) 4l|: 



cp{t) ^ I d'x{p^{x)p'^m = 



cf^^ 



F2 S 



E2 4 



2S2 ASPQ(/i,,/i,) 



F2 



1 



2 S 
S2 4 



Th{t/T), 



(8) 



- 



Th{t/T). 



(9) 



The operators are separated in time by t and projected onto zero spatial momentum. At 
the NLO an "effective" chiral condensate Sofr, 

3/5i 



^efr 



S 1 + 



appears with a combination ^f^ = niiJ^esV. The parameter /3i is the so-called shape coef- 
ficient that depends on the shape of the box. In our numerical study, /3i = 0.0836. In Eqs. 
dH]) and dnD, the time dependence is simply written by a quadratic function hi(t/T) as 



hit/T) 



t 

f 



1 

12 



(11) 



Next, we consider the flavored axial-vector and vector operators Aq{x) = q'(a;)r"7o75g(x) 
and Vq{x) = g(x)r'^7og(x). These correlators to O(e^) in the Nf = 2 unquenched ChPT are 



39| 



CA{t) ^ j d'x{Al{x)Aim = I J° + |j {y^.Jl - yk,oJ^ 

4/i,sr(/x.)T2 



F2 



V 



hiit/T) , 



(12) 



Cv{t) 



d'x{V,^{x)V,^m = I J° + A (^_^ JO _ ^k,oJ^^ } , (13) 

where koo is another numerical factor depending on the shape of the box. In our numerical 
study, koo = 0.08331. jT^'s are defined by 



= 3 I 2 + 2 



S 



J' 



4-2 



sl-"(/.f)V , 1 sr'(/xf) 



^fuU^^eff-j \ 2 -I vfull/',,efr 



cff 



(/^ 



, (14) 
■ (15) 



It should be noted that the axial-vector and vector correlators are sensitive to F, as it 
appears as an overall constant. Their time dependence and other factors represent the NLO 
effect. For the pseudo-scalar and scalar correlators, on the other hand, S determines the 
overall constant and F appears only in the NLO correction term. Therefore, by using both 
types of correlators the two LECs can be extracted with a good sensitivity. 



III. LATTICE SIMULATIONS 



We summarize the setup of our numerical simulations. Details of the configuration gen- 
eration and the eigenvalue calculations are given in 



Our lattice is 16^ x 32 at a lattice spacing a = 0.1111(24) fm determined from the Sommer 



scale tq = 0.49 fm as an input 
Dirac operator 

Dim) - 



321 ]. We employ the overlap fermion 28|, |29|], defined by the 



mo + 



m 



m 



mo 



75Sgn[ifiy(-mo)] 



(16) 



for a quark mass m. Hy/{—mQ) = '-f^Dwi—mo) denotes the standard Hermitian Wilson- 
Dirac operator at a large negative mass —mo. We choose mo = 1.6 throughout this work. 
We remark that the mass parameters are given in the lattice unit unless otherwise stated. 



We employ the Iwasaki action |33l . 



34| for the gauge field at /? = 2.35 with an additional 



determinant factor det[Hyy/[H^r + m f)] i n the partition function produced by extra Wilson 
fermions and twisted-mass ghosts SS,^, Is?] (we set mt = 0.2). With this choice the global 
topological charge does not change its value during the molecular dynamics updates of the 
Hybrid Monte Carlo algorithm. Fixing topology in this way is desirable for the study of the 
e-regime of QCD, since the analytical expressions of ChPT are given at a fixed topological 
sector. In this work, we take the trivial topological sector u = only. 

The sign function in ( fT6i) is approximated by a rational function with Zolotarev's co- 
efficients after projecting out a few lowest eigenmode's contribution. With 10 poles the 
accuracy of the sign function is lO"*^^"'^-'. Therefore, the violation of the chiral symmetry 
due to the lattice action is negligible in our work. 

In this work, we use the gauge configurations generated at sea quark mass m = 0.002, 
which corresponds to ~ 3 MeV. For this value the parameter mHV is about 0.556, and 
the system is well within the e-regime. As we decrease the quark mass from the p-regime 
(mHV ^ 1) to the e-regime {mT,V < 1), the lowest eigenvalue of D{my D{m), ~ + m^. 
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is bounded from below by a repulsion of the lowest eigenvalue from zero due to the fermion 
determinant Yiki-^l'^''^'^)- The condition number of the operator D{myD{m), and thus the 
computational cost for its inversion, saturates near the boarder between the p-regime and 
the e-regime. As a consequence, for example, the number of multiplication of the Wilson- 
Dirac operator needed per trajectory at m = 0.002 is only about 1.5 times greater than 
that at 10 times heavier sea quark mass 18|. We have accumulated 4,600 HMC trajectories 
after discarding 400 trajectories for thermalization. The numerical cost is about one hour per 
trajectory on a half rack (512 nodes) of the IBM BlueGene/L (2.8 TFlops peak performance). 

At every 10 trajectories, we calculate the meson correlators at four values of valence 
quark masses m = 0.0005, 0.0010, 0.0020, and 0.0030 ranging 1-4 MeV in the physical unit. 
For the inversion of the overlap-Dirac operator we use the multi-shift Conjugate Gradient 
(CG) solver to calculate all the valence quark propagators simultaneously. The solver is 
accelerated by projecting out the subspace spanned by 50 pairs of lowest-lying eigenmodes. 
(The eigenvalues of -D(O) form a pair with their complex conjugate; the eigenvector of the 
counterpart is produced by multiplying 75.) With this projection, the solver performance is 
an order of magnitude better and roughly independent of the quark mass. These eigenmodes 



of -D(O) has been calculated using the implicit 
disks for studying the eigenvalue distribution 



y restarted Lanczos algorithm and stored on 



17 



18|. 



We compute the meson correlators with the Low Mode Averaging (LMA) technique 
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451 ]. Using the Nep pairs of lowest-lying eigenmodes, we decompose the quark propagator 
D(m)~^ into the low-mode contribution [-D(m)];~^ and the rest [D{m)]l^lg,^ as 



Dim)-\x,y) = [Dim)]-l{x,y) + [D{m)]-l^,{x,y) 

^ Ukix)ul{0) 75Mfc(a;)Ml(0)75 



E 

k=l 



+ 



m/2mo)Xk + rn {1 — m/2mo)Xl + m 



[DM]-l^,ix,y), 
(17) 



where Uk{x) is the eigenvector of -D(O) associated with its eigenvalue A^. While the high- 
mode contribution have to be obtained by using the CG solver for a fixed source point y, 
the low-mode contribution can be calculated from the 50-1-50 low-lying modes for any source 
and sink points without extra cost. With the LMA technique, we average over the source 
point for a part of the meson correlator that is purely composed of [D(m)]j'J^. For other 
contributions, we simply use a fixed source y at the origin. 
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FIG. 1: Pseudo-scalar (left) and scalar (right) correlators with (solid curve) and without (dotted 
curve) the low-mode averaging. Correlators are calculated on a single gauge configuration at 
m = 0.002. 

The LMA technique is effective to improve the statistical signal when the correlator 
of interest is dominated by the low-mode contribution. In Figure [T] we demonstrate the 
improvement by taking pseudo-scalar and scalar correlators calculated on a single gauge 
configuration as an example. Some wiggle observed without LMA is completely washed out 
with LMA, and smooth curve is obtained. A similar comparison is shown in Figure [2] for the 
axial-vector and vector correlators, for which the improvement with LMA is marginal. This 
indicates that these correlators are not simply dominated by the low-lying modes. Note that 
the magnitude of these correlators is two orders of magnitude smaller than pseudo-scalar 
and scalar correlators. 

For the statistical analysis, we use the jackknife method with a bin size 20, which cor- 



error saturates for the 
. Since the low-lying 



responds to 200 HMC trajectories. With this choice the statistical 
calculation of the averaged low-lying eigenvalues (A^) as studied in 
modes reflect the long-distance physics (zero- modes of pion fields in the language of ChPT), 
these quantities are expected to have the longest auto-correlation time among other physical 
quantities. In Figure [3l we plot the Monte Carlo history of the quantities of interest, i.e. 
the pseudo-scalar and axial-vector correlators at the largest time separation t = T/2 = 16, 
Cp{T/2) and Ca(T/2) (see the definition given later). We observe that the 200-trajectory 
gives a reasonable range of the auto-correlation. However, the integrated auto-correlation 



time calculated following the definition in 46| is substantially shorter: 37(12) and 47(16) 
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FIG. 2: Axial-vector (left) and vector (right) correlators with (solid curve) and without (dotted 
curve) the low-mode averaging. Correlators are calculated on a single gauge configuration at 
m = 0.002. 
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FIG. 3: The Monte Carlo history of the pseudo scalar correlator (top) and the axial vector 
correlator (bottom) at t = T/2 = 16. 



HMC trajectories for Cp{T/2) and Ca(T/2), respectively. We therefore conclude that the 
binsize of 20 is a conservative choice. 
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FIG. 4: Axial-vector correlator in the e-regime. Filled square symbols denote the low-mode 
averaged correlator while the crosses are not averaged. The solid curve shows the best fit with the 
ChPT formula. 

IV. NUMERICAL RESULTS 



In this section we explain the fit of our data to the (partially quenched) ChPT formulae 
in Eqs.®, ©, and 

First, we study the axial- vector current correlator at = = 0.002 (full QCD point). 
Since we use the local axial current Aq{x), which is not a conserved current on the lattice, we 
need a finite renormalization to relate the lattice current to the continuum current ^o(^) 
^o(^) = ZaAq{x). We calculated the renormalization factor Za non-perturbatively through 
the axial Ward Identity and obtain Za = 1.3513(13). In the following, the values quoted for 
F include this Z-factor, neglecting its tiny statistical error. 

Figure m shows the axial- vector current correlator with and without the LMA technique. 
Although the improvement by LMA is marginal for this channel, the statistical error is re- 
duced by about 30%. We fit the data to the ChPT formula (fT2l) taking F and S as free 
parameters. With a fitting range t G [12,20] we obtain the solid curve shown in Figure HI 
Roughly speaking, the overall magnitude (the constant piece) determines F, while the cur- 
vature (or the term proportional to hi{t/T)) gives S. We obtain S ~ [260(32) MeV]'^ and 
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FIG. 5: Vector correlator in the e-regime. The sohd curve is the ChPT prediction with inputs 
S = [227.6(3.7) MeVf and F = 87.3(5.6) MeV (No free parameter left). 

F ~ 90(6) MeV. The statistical error is large for S (about 30%), because it is extracted 
from a tiny curvature. 

A similar plot is obtained for the vector channel as shown in Figure [51 At the leading 
order, i.e. in the massless limit at finite V, axial- vector and vector correlators become 
identical because of the exact chiral symmetry; the difference arises due to the zero-mode 
integrals when /i is finite. The line in the plot shows the formula (|T3|) with the parameters F 
and S obtained from the fit of the axial- vector and pseudo-scalar channels (see Section IVT). 
It shows a remarkable consistency. 

Next, we consider the pseudo-scalar channel at m„ = = 0.002. The lattice data ob- 
tained with the LMA technique are shown in Figured With F obtained via the axial vector 
correlator as an input, we compare the pseudoscalar channel with the chiral perturbation 
theory formula at the NLO ([8]). The motivation for this choice is that the pseudo-scalar 
channel itself does not have a good sensitivity on F, as it appears only at the NLO. From 
a fit in the range t G [12, 20] we obtain the chiral condensate as S ~ [227(4) MeV]^. 

A similar result can be obtained for the scalar channel as plotted in Figure [71 We draw a 
curve representing the ChPT formula with the parameters obtained from the fit of the 
axial- vector and pseudo-scalar channels (see Section [V]). Again for this channel, the data 
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FIG. 6: Pseudo-scalar correlator in the e-regime. The solid curve represents a fit with the NLO 
ChPT formula. 
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FIG. 7: Scalar correlator in the e-regime. The solid curve represents a fit with the NLO ChPT 
formula. 
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are consistent with the expectation within the statistical error. 

Using the data at 7^ we can check the consistency of the lattice results with the 
partially quenched ChPT formulae ([S]) and (jH]). Figure [H] shows the pseudo-scalar and scalar 
correlators at four different valence quark masses ranging 1-4 MeV while fixing the sea quark 
mass at 3 MeV. The shape of the correlators does not strongly depend on the valence quark 
mass, but some dependence can be seen. With the input parameters F and S determined 
from a global fit of the axial- vector and pseudo-scalar channel (see Section |V]), we draw 
the expectation from the partially quenched ChPT in Figure [H There is no additional 
parameters in this analysis, hence it gives a stringent test of either the NLO ChPT formula 
or the lattice calculation. We find an excellent agreement for all the valence quark masses 
within the statistical error. Some specific points of the correlator, e.g. at t = 12, 14, or 
16, are plotted as a function of the valence quark mass in Figure [H The decrease of the 
correlators with the valence quark mass is nicely reproduced by the lattice data. 

V. EXTRACTION OF THE LOW ENERGY CONSTANTS 

Since the axial-vector and pseudo-scalar channels are complementary to each other in 
determining the LECs F and S, we fit both channels simultaneously with these two param- 
eters to obtain our best results: S = [227.6(3.7) MeV]^ and F = 87.3(5.6) MeV. Here, the 
statistical error of the lattice spacing a is also taken into account. The fit range is again 
t G [12, 20] for both channels. Under the change of the lower limit of the fitting range from 
10 to 15, which corresponds to 1.1 fm to 1.7 fm, the fit results are quite stable (within 1%) 
with similar error bars. 

Multiplying the non-perturbative renormalization factor to convert our result of the 
chiral condensate to the continuum MS scheme we obtain E'^^(2 GeV) = [239.8(4.0) MeV]^, 
where the error represents the statistical one. 

Possible sources of the systematic error include the discretization effect of O(a^) and the 
higher order effect of the e-expansion. Since we do not have data at different lattice spacings, 
we cannot quantify the discretization effect beyond a rough order counting. Assuming that 
the relevant physical scale is the QCD scale Aqcd ~ 300-500 MeV, the size of (9((aAQCD)^) 
effect is 3-8%. This small scaling violation is supported for the overlap fermions by recent 
quenched simulations (in the p-regime) albeit for different physical quantities. The 
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FIG. 8: Partially quenched pseudo-scalar (top) and scalar (bottom) correlators in the e-regime. 
The solid curves are the ChPT results with S = [227.6(3.7)MeV]2 and F = 87.3(5.6)MeV as inputs 
(No free parameter left). 

unphysical (heavy) Wilson fermions we introduced to fix the topological charge should not 
have any negative impact on the scaling, because they never arise in the external states and 
merely affect the gluon action at order a? in the Symanzik's effective theory. Their effect on 
the effective gluon action is further minimized by the accompanied ghosts that cancel the 
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FIG. 9: Valence quark mass dependence of the partially quenched pseudo-scalar (top) and scalar 
(bottom) correlators in the e-regime at each time slice. 

Wilson fermion contributions except for the near-zero modes of H^, that appear for locally 
bumpy gauge configurations, or for the so-called dislocations js^. 

The higher order effect of the e-expansion appears due to the finite volume lattice. In our 
analysis, the NLO terms, i.e. 0{e^) terms, are included and the remaining corrections are 
of 0{e'^) ~ C(l/(Acut-^^)'^), with Acut the cutoff scale of ChPT. With a conservative choice 
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Acut ~ 300-500 MeV, this uncertainty is 0.3-2%. This order-counting gives a reasonable 
estimate for the NLO correction to S. Namely, the size of the NLO correction (See — S)/E 
(see f llOp ) estimated with this order-counting is 14%, while the real correction was 20%. This 
agreement indicates that the order-counting argument is indeed valid up to an 0(1) factor, 
which depends on the shape of the space-time box 6| and could be sizable at the NNLO. 
At the same order of the e-expansion, there is a contribution from the finite pion mass 
(9((m7r/Acut)^), which numerically gives about 2%. We note that such a small uncertainty 
does not apply for the curvature of the correlators, since the curvature itself is the quantity 
arising at the NLO. But the extraction of F and S relies mainly on the constant piece of 
the correlators, that starts from the leading order in the e-expansion. 

Adding the uncertainties from the scaling violation and from the higher order effects in 
the e-expansion in quadrature we estimate the dominant systematic error is of order 9%. For 
more robust estimate beyond the order-counting, it is important to study the finite lattice 
spacing and volume effects in the future works. 

The result for the decay constant F may be compared with the calculation in the p- 
regime Ml]- In 52] a preliminary result for the pion decay constant with two-flavors of 
dynamical overlap fermion is presented. With the same lattice action, the calculation is 
done at slightly coarser lattice spacing, a ~ 0.12 fm. An analysis with the next-to-next-to- 
leading order (NNLO) ChPT yields F = 78.6(2.7) MeV (the error denotes the statistical 
one), which is slightly lower than the calculation in this work. Assuming the 9% systematic 
error, however, both results are consistent with each other. 

The result for the chiral condensate S may be compared with the extraction from the 
lowest-lying eigenvalue through the chiral Random Matrix Theory HQ: S^^^(2 GeV) = 
[251 (7) (11) MeV]^. Here, the errors represent the statistical and the systematic due to the 
NLO correction in the e-expansion. The calculation in this work from the meson correlators 
includes the NLO correction, and in fact the result deviates from the previous result by about 
the estimated NLO correction. When one integrates out the non-zero modes of the pion field 
from the chiral Lagrangian, there remains the zero-mode integrals with an "effective" chiral 
condensate, which is given as 

at the one-loop order. The one-loop correction is a substantial effect 
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53j: a factor of 1.202. 



If we correct the value in this work with this amount we obtain [255(4) MeV]^, which is now 
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calculation has been given recently 
the topg 
(See 



1 statistical error. This remarkable consis- 



24| , and an argument based on an analytical 
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5J] . It is also notable that our results are consistent with 
ogical susceptibility, from which we extracted S^^(2 GeV) = [254(5)(10) MeV]^ 



56| for the details). 



VI. CONCLUSION 



The extraction of the low-energy constants F and S from the meson correlators in the 
e-regime has a significant advantage over the conventional approach. Already at the NLO 
in the e-expansion, the remaining higher order effect is a per cent level, and thus a precise 
calculation of both F and S is possible without delicate chiral extrapolations. 

The dynamical lattice simulation is feasible with a small sea quark mass ~ 3 MeV on 
a 16'^ X 32 lattice with L ~ 1.7 fm, where the scale is determined assuming tq = 0.49 
fm. Meson correlators are obtained with a good precision by using the low-mode averaging 
technique. The numerical results passed all the consistency checks for different channels 
as well as for the partially quenched mass combinations. Using the NLO ChPT formulae 
in the e-expansion, we obtain S^(2 GeV) = [239.8(4.0) MeV]^ and F = 87.3(5.6) MeV 
with a possible systematic error of order 9% which is dominated by the discretization error. 
Reducing the systematic error by repeating the calculation on a larger lattice at smaller 
lattice spacing will be an important future work. 
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